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We provide the Barrett-Crane spin foam model for quantum gravity with a discrete action princi- 
ple, consisting in the usual BF term with discretized simplicity constraints which in the continuum 
On . turn topological BF theory into gravity. The setting is the same as usually considered in the lit- 

■ erature: space-time is cut into 4-simplices, the connection describes how to glue these 4-simplices 

' together and the action is a sum of terms depending on the holonomies around each triangle. We 

impose the discretized simplicity constraints on disjoint tetrahedra and we show how the Lagrange 
5_( . multipliers distort the parallel transport and the correlations between neighbouring simplices. We 

then construct the discretized BF action using a non-commutative ★-product between SU(2) plane 
\ waves. We show how this naturally leads to the Barrett-Crane model. This clears up the geometrical 

. meaning of the model. We discuss the natural generalization of this action principle and the spin 

T-H ' foam models it leads to. We show how the recently introduced spinfoam fusion coefficients emerge 

CO ' with a non-trivial measure. In particular, we recover the Engle-Pereira-Rovelli spinfoam model by 

, weakening the discretized simplicity constraints. Finally, we identify the two sectors of Plebanski's 

' ^ ' ' theory and we give the analog of the Barrett-Crane model in the non-geometric sector. 

cr: 

Introduction 

m ■ 

^ ] Striving to make sense of a quantum theory of gravitation, several approaches have converged towards the con- 
struction of the spin foam formalism as a promising non-perturbative framework to describe quantum theories in a 
background-independent context. A spin foam is ultimately a combinatorial object given by a 2-complex, possibly 
^■f-^ ■ seen as dual to a triangulation, whose simplices are labeled by data from the representation theory of a given group. 
If loop quantum gravity has given a clear picture of the kinematics of the theory in terms of spin networks, spin foams 
are expected to provide a representation of their (gauge) dynamics, generated by the Hamiltonian constraint, and 
' thus to give a tool to define transition amplitudes [1| . On the other hand, it has been shown that path integrals for 
, lattice gauge theories can be formulated as sums over spin foams living on the lattice In particular, spin foams 
living on a single triangulation are well adapted to topological BF theories P, 0I- In the 3d case, where SU(2) BF 
theory corresponds to gravity, they have been shown to yield the correct quantization The group representation 
data defining each spin foam can be interpreted, using the canonical analysis, as quantized flux observables @. 
5— i Several models have applied for describing quantum gravity. Among them, the most studied is the Barrett-Crane 

^ (BC) model 0, the first non-topological model. The study of its semi-classical limit has revealed that the two-point 
function exhibits the correct behaviour with regards to the distance Q, but that it fails to reproduce the expected 
tensorial structure of the free graviton propagator of spin 2 [§]. New propositions have been made to correct the 
drawbacks of the BC model [lOl . [Til . [T^ . \\^ , including a Lorentzian version |14| | , and with a finite Immirzi parameter 
[Tst - The edification of these works is based on different foundations. First, there exists a reformulation of the action 
of general relativity as a topological BF theory with additional constraints C quadratic in the 2-form B G so (4): 



X 



Sgr{B,A)^Sbf{B,A) + \C{B) 



(1) 



where A stands for the Lagrange multipliers imposing these so-called simplicity constraints. They express that B 
is the wedge product of a tetrad, restoring the degrees of freedom of GR flGj. The idea is thus to start with the 
standard spin foam model for BF theory and to introduce the constraints at the quantum level. Their discretization 
involves different subsimplices of the triangulated manifold, the most discussed being the so-called diagonal and cross- 
simplicity constraints, respectively taking place at each triangle and each tetrahedron. Efforts have thus concentrated 
on the quantization of the geometry of a tetrahedron , which is the natural arena for the Hamiltonian framework, 
using geometric quantization [l^ to define a state space and quantum observables. According to the syrnplectic 
structure, the variables discretizing the field B are promoted to generators of so (4) at the quantum level (see [ill [Tsj 
typically). The constraints can then be solved inducing restrictions on the spin data labelling the quantum states of 
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the tetrahedron. From this Hamiltonian perspective, the BC model strongly imposes the constraints. But they form 
a second class system (see [l9| for a comparison of the different models in this regard) , and have thus to be weakly 
imposed to avoid the loss of degrees of freedom. This is precisely the purpose of [I^l ■ 

The issues are then to insert these results into the spin foam and to stick tetrahedra together. This is usually 
done by identifying the spin data of the spin foam with the labels of the tetrahedron state space, as in Q in the 3d 
case. Indeed, it enables to assign the same spin to each copy of a triangle as a regluing process. Notice also that 
the amplitude of BC model was built considering a single 4-simplex. From the point of view of [l^ . the regluing of 
tetrahedra is correctly performed in this model, but that of 4-simplices is failing, forgetting that each tetrahedron is 
shared by two 4-simplices. 

The standard way thus involves a recourse to both the Hamiltonian and the Lagrangian methods, together with 
a gluing process which takes place at their meeting point. We here consider the purely Lagrangian point of view, 
discretizing the functional integral. This approach has already been explored to perturbatively define and compute 
generating functionals of theories written as a BF term plus a (local) polynomial potential in B [2d[ . This work is in 
particular based on the ability of performing the integration over the field B in pure BF lattice model. It remains 
unclear how such an integration can be exactly performed for a constrained BF theory like ([T]) (see however [2l| 
for an interesting path integral formulation of the new spin foam models, but not exactly based on the standard 
BF discretized action). A second aspect in [2^] and i^lj, which is certainly related to the ability of performing the 
integration over B, is the breaking up of the simplicial decomposition into disjoint 4-simplices. One thus consider 
an action for each 4-simplex and compute its partition function in terms of boundary variables for the connection, 
without an explicit regluing for B (except the above-mentioned process based on the canonical analysis). 

In this note, we provide an action for the BC model, which consists of a BF term supplemented with the usual 
discretization of the simplicity constraints. This makes clear the physical setting of the BC model. The BF part of 
the action is built on disjoint tetrahedra, and the constraints are independently imposed on each of them. There is 
no explicit regluing of the variables B assigned to different copies of a triangle, but a gluing process with boundary 
variables for the connection instead. Such a process is obviously sufficient for pure BF theory. However, we will 
show that the Lagrange multipliers imposing the contraints behave as sources for the curvature, and make this gluing 
fail to give the correct relations of parallel transport. The BC model has long been suspected because of its lack of 
correlation between neighbouring simplices. It clearly appears here that the key point in this regard is the gluing 
between tetrahedra once simplicity has been taken into account for each of them. As a way-out in the spirit of 
Regge calculus, one may think of a gluing involving the field B, before its integration. We also provide a very simple 
parametrization of these correlations which could be useful to get some flexibility to modify the BC model. 

We clarify some issues concerning the BC model. Our parametrization of the constraints enables to identify the 
physically relevant part of the partition function and to use a freedom in the measure to show that the BC model 
weakly implements simplicity with specific fluctuations. The measure to be used is the most natural from the point 
of view of SU(2) lattice gauge theory. Another essential aspect is that the formulation ([1]) contains gravity and also 
a non-geometric sector. The two sectors are clearly identified and a spin foam model for the analog of the BC model 
in the non-geometric sector is g iven. 



The spin foam models of [iCt |1 2|| intend to correct the BC model in the view of different defects: it is supposed 
to strongly impose simplicity for~|lo|. or to be built on disjoint 4-simplices for p^ . The main feature of these models 
is the amplitude of a 4-simplex, based on the 15j-symbol for Spin(4) with some fusion coefficients which supplant the 
BC intertwiner. The amplitudes for triangles and tetrahedra however differ. The fusion coefhcients are here shown 
to appear in the framework of the BC model with a non-trivial measure. They are thus very natural objects when 
trying to implement metricity and are not specific to the new models. 

The organization of the paper is as follows. The first section reviews standard results of spin foam quantization 
for Spin(4) BF theory. In the second section, we study the implementation of diagonal simplicity in the functional 
integral. This is indeed usually considered as the starting point to impose cross-simplicity, though it has never been 
precisely studied from this point of view. This is the simplest modification of the BF spin foam model making it 
non-topological. We tackle the issue of cross-simplicity in section IIIII The analysis is based on the formulation of 
P^, and the BC model is derived from a discrete action principle. We use a freedom in the measure to generalize the 
BC model in section HVl and show how the fusion coefficients of the new spin foam models naturally emerge in this 
setting. The measure is also used to weaken the implementation of the constraints so as to recover the new spin foam 
models. We close this note with a discussion concerning the physical setting of the BC model and its improvement. 

During the process of writing up the present paper, Conrady and Freidel [2ll . [22| also provided a path integral 
formalism for the new Engle-Pereira-Rovelli and Freidel-Krasnov spin foam models. Their approach is different from 
ours since the purpose of their path integral is to exactly reproduce the amplitudes of these spin foam amplitudes, 
while our formalism starts from the generic path integral for the discretized BF theory and explores how to impose the 
simplicity constraints in a clear geometric way. At the end of their day, the present work allows to derive a larger class 
of spinfoam models and provides them with a clear geometric meaning while the Conrady-Freidel framework is more 
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specific to the Engle-Pereira-Rovehi and Freidel-Krasnov spin foam models and allows to extract their (semi-)classical 
behavior. 



I. REVIEW OF SPIN FOAMS FOR BF THEORY AND THE DISCRETIZED SIMPLICITY 

CONSTRAINTS 

Let us first review some basics facts about the conventional spin foam model for topological Spin(4) BF theory. 
We will all along use the sphttings of the group, Spin(4)= SU(2) x SU(2) and of its algebra, spin(4) — su(2) © su(2). 
Consider a triangulation of a 4-manifold. The basis of Regge discretisation, in four dimensional spacetime, is to 
concentrate curvature on 2-simplices. Tetrahedra and 4-simplices are thus flat, and immmersed in the 4-manifold 
with arbitrary, and so a priori different frames, that is local coordinates for the 'internal' space, carrying the vectorial 
representation of Spin(4). The spin foams then lie on the 2-complex dual to the triangulation. We will indifferently 
denote the triangles and their dual faces by /, the tetrahedra and the edges of the dual faces by r, and the 4-simplices, 
dual to vertices, by v. 

In the continuum, the theory is defined by the classical action Sbf = J ti{B A F), where i? is a 5pin(4)-valued 
2-form, and F the curvature of a Spin(4) connection. We define the generators = —i(Ji/2 of 5u(2), where the 
matrices ai are the Pauli matrices, which are anti-hermitian matrices satisfying [Ti,rj] = e^j r^.. The trace 'tr' over 

su(2) is the conventional trace, with the normalisation: tr(TiTj) = — The discrete connection consists of Spin(4) 
group elements Gyr — (g+y-n g-vr) allowing parallel transport between frames along the boundary edges of the dual 
faces. The orientation of each dual edge induces an orientation between each pair of 4-simpliccs sharing a tetrahedron. 
The basic variables are, for each r, the holonomies Gs(t)t a-nd Gt^t) representing the parallel transport between the 
4-simplex source s(r) and r, and between r and its 4-simplex target ^(t). Composing these elements and given 
an orientation of each face, we define the holonomies around the faces^, Gf{vQ) = GygTiGnviGy^r^ ' ■ ■ > ^^o being a 
base-point for /. Such a discrete holonomy can be transported in another frame of the boundary of /, e.g. that of a 
tetrahedron, by the following obvious formula: 

Gf{T)^G-^,Gf{v^)Gy,r (2) 

where Gy^r turns around / using the onwards path from vq to r. We discretize the continuous 2-form B on the 
triangles of the triangulation, with respect to the local frame of a tetrahedron, Bf{T), or that of a 4-simplex, Bf{v). 
Each Bf is a, bivector, that is an element of spin(4), Bf = 6+/ ® where each b±f can be seen both as a su(2) 
matrix or a 3- vector in agreement with its definition: b± — V±Ti, with Vj^^ — \{\^^ jk^''^ i B'^^). There is several 
ways to consistently define the bivectors together with the rules for parallelly transporting them. A natural way^is 
to first define Bf{vo) for a base-point vq, dual to a 4-simplex, and then, using the orientation of / and the elements 
Grt), define the bivector in the other frames: 

Bf{ri) = G-^\ Bfivo) Gy„r, , Bf{v,) = G;^ Bf{n) Gr,y, (3) 

and so on until t„ is reached as in figure [TJ The relations ^ and ([3]) ensure that the discrete action for BF theory, 
shown below, does not depend on the choice of a base-point for each dual face. 

The Spin(4) gauge covariance reflects the arbitrariness in the choice of local orthonormal bases for the (flat) 
4-simplices and tetrahedra. The gauge transformation K thus acts through Spin(4) group elements Ky and Kt 
associated to 4-simplices and tetrahedra which rotate the local frames, as in the framework introduced in [loj : 

K > Bf{v) = Ky Bf{v) (4) 

Kt>Gyr^ Ky Gyr K^^ (5) 

The gauge invariant discretized BF action is then: 

= ^Tr(i?/(z;o)G/(«o)) =Et'^(^+/(«o)-9+/(«o)) (6) 



^ We suppose for the sake of notation that the orientation of each dual edge coincides with that of the dual face /. 

^ Another method often used in spin foams is to divide dual faces into several pieces corresponding to the choice of a boundary simplex. 
The so-called wedges correspond to pairs {f,v). Geometrically, it means that the initial triangulation has been broken up into disjoint 
simplices. One can then define independent bivectors for each piece and reglue simplices in a specific way. It will be detailed in the next 
sections, for we will argue that the natural setting of the BC model corresponds to using pairs (/, r). 
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Bf{vo) 




FIG. 1: The dual face / is equipped with an orientation and a base-point. The bivector _B/(wo) is defined in the frame of this 
base-point, and then parallelly transported along the boundary of / until Bf{T„) is reached. This ensures that the action for 
BF theory does not depend on the choice of the base-points. 



in which the traces are independent of the chosen frames thanks to It is also clear that ^ is invariant under 
the change -B/(wo) = C!f{vny^Bf{vo)Gf{vo)~" for any n, since: ti{B f{vo)G f{vo)) = tr{Bf{vo)Gf{vo))- This way of 
thinking avoids imposing that Bf(vo) is left invariant by G/(wo). Such a relation, well-known in Regge calculus, comes 
in fact from the discretization of an equation of motion of Plebanski's action (we come back to this point in sections 
nil CI and lV|) . Thus it may be sufficient to get it as e.o.m. for the discrete action principle as in for the new spin 
foam models. 

The spin foam model is then defined as the partition function for this action, and can be interpreted as a discrete 
way to compute the volume of the moduli space of flat connections. Indeed, the dynamical variables are the bivectors 
and the holonomies along the dual edges Gr = Gst^T^T-GT-ti^T--^, and the partition function reads: 

Zbf = j\{dGr j \{dBs e'i:/Tr(G,HB,w) (7) 

r ■> f 

\{dGr \{6{Gf) (8) 

/ 

where dB is the Lebesgue measure on spin(4) = ® R'^, and dG — dg+dg- is the Haar measure on the group 
Spin(4)= SU(2) x SU(2). The action being linear in each _B/, their integrations lead to delta functions over the 
group for the holonomies around the dual faces'^ . This is precisely the types of computations we want to perform for 
non-topological theories. In particular, the result of the integrations over the bivectors expressed as a function over 
the group gives a great part of the physical content of the theory. Here, it is a simple discretization of the condition 
asking for flatness of the connection [3] and, as in the continuum, bivectors are straightforwardly interpreted as 
Lagrange multipliers. Expanding the condition S{Gf) according to the usual formula 6{g) — djXjig) for SU(2) 
group elements g, where dj = 2j + 1 is the dimension of the irreducible representation of spin j, the dual faces are 
labelled with representations , j_/): 

Zbf= Y[dGr U^j+fdj^^ Xj+f{9+f)Xj-f{9-f) (9) 

In the continuum, the action of general relativity can be written as that of BF theory supplemented with additional 
constraints which ensure that the field B comes from the wedge product of a tetrad. They read: 

eijKL B^^l B^^ cx e^^Acr (10) 



^ The integrals over Bf is not altered when changing f or g_f into its opposite, so that the delta functions are in fact over SO (3) X SO (3) 
(see next section). 
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where fi, v, X, a are spacetime mdices. In fact, the theory consists m two sectors, only one describmg gravitation, 
which we wih cah the geometric sector, with B^'^ = e^'^ j^i^ A for a tetrad 1-form e. The non-geometric sector, 
often given the misleading name of topological, is characterized by B^"^ = A e'^ (see [l^ for details). Its equations 
of motion only impose the vanishing of torsion and leave curvature free. The so-called simplicity constraints are 
conventionnally discretized on simplices of different dimensions: 



i)Diagonal simplicity constraints: cuklB ^ B ^ — 0, 
ii)Cross-simplicity constraints; tuKhSy B^,^ — 0, 
iM)Volume constraints: eijxLBy B^,^ cx T4 



for a single triangle (11) 
for two triangles of a tetrahedron (12) 
for / and /' sharing a point in the 4-simplex v (13) 



in which we have dropped the frame labels. The diagonal constraints express that the bivector of each face is simple, 
while the cross-simplicity constraints relate the bivectors of each tetrahedron. The third constraint (fT5)) , taking place 
at the level of each 4-simplex v, ask for the volume of v to be independent of the bivectors used to compute it. It 
has been shown (see fld\ for instance) that can be replaced with the closure relation, X]/car ^/(''') ~ ^ each 
tetrahedron. The latter is usually taken for granted since it is an equation of motion of BF theory. It will also be the 
case here for the models under consideration. 

Thinking of adding the constraints to the BF action ([6]) , a natural idea is to use the fact that they are quadratic in 
the field B, so that Gaussian integrations can be performed. However, this does not turn out to be a great advantage 
in this context, even for the diagonal constraint pT|) which involves only one triangle. Indeed, expressions pTjl . IT^ 
and (jl3p are not usual polynomial integrations in B, but constraints imposed with Lagrange multipliers instead. One 
has then to deal with determinants involving these Lagrange multipliers. This approach can be nevertheless fruitful 
in other contexts, for example in [27] which investigates the renormalization of Plebanski's action. A reformulation 
of the discretized constraints PT|) and has been introduced in llOj] , and used in [ll|, [l^, [11] , which clarifies their 
geometric interpretation. It simply expresses that a tetrahedron whose triangles are labelled by bivectors satisfying 
the constraints spans a three dimensional subspace of Minkowski space (that is the 'internal space'), so that there 
exists for each tetrahedron a unit vector Nr orthogonal to each of its faces in the following sense: 

eijKL N:^Bf'^ ^0 or N^jBy^O, (14) 

respectively for the gravitational and non-geometric sectors. In the following sections, we will use another parametri- 
sation, using SU(2) group elements, which can be straightforwardly interpreted in terms of ([H]). 



II. SIMPLICITY OF A SINGLE BIVECTOR 



In this section, we introduce basic ingredients from which non-topological spin foam models can be derived in a 
clean way, by looking in details at the diagonal constraints pip . First, it can be seen as an exercise to the challenge 
of imposing the cross-simplicity constraints. Our derivation of the BC model will be based on the same method. 
Second, the effects of the diagonal constraints have not been investigated so far from the point of view of the lattice 
path integral. The physical content of pT|) is that each Bf, or equivalently their Hodge dual, *By = xL^f^ ^ 
is simple, that is to say that it is an anti-symmetrized product of vectors, u'^f '^^ . It thus corresponds to a genuine 
non-topological model. 

The diagonal constraints pip are usually imposed thanks to a geometric interpretation which identifies the Spin(4) 
representations (j+/j J-/) labelling faces in the state-sum model ^ for BF theory, with a quantization of the triangle 
areas. Such an identification is motivated by the canonical analysis of the discrete action. It then leads to implementing 
the diagonal constraints by imposing the equality of the self-dual and anti-self-dual representations j+/ = 
which corresponds to the simple representations of Spin(4). Here, we show in a rigourous way to what extent this 
implementation is correct. 

The usual imposition of the constraints pip , which borrows tools from geometric quantization, leads to the following 
partition function: 

^simple = jWdGr WZiiGf), (15) 

where Zf{Gf) = ^ d% x,, (ff+z) X„ (s-/) (16) 

with k = 2 defining the usual face amplitude. We will show how the discrete path integral for the BF action with the 
diagonal constraints naturally gives p6p with fc = 0, that is a trivial face amplitude, and how this result generalizes 
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to: 

ZfiGf) = x{9+f) ^if (d+f) (s-f) (1"^) 

where x is a SU(2) class function, typically a character, to ensure the gauge invariance of Zf. x is a measure ambiguity 
of the path integral, and as we will see, the models defined by (flT)) share the same basic physical features due to (|lip . 
When X is a character, it induces simple shifts of the spin of the self-dual character. The reason why ITTl) does not 
break the symmetry between the self-dual and anti-self-dual sectors will be clear below. 

Since our aim is to explicitly perform the integrals over the bivectors in the partition function, the key point is to 
find a convenient parametrization of the solutions to (|lip . Decompose the bivectors Bf into self-dual and anti-self-dual 
parts: Bf = 6+/ ® ^-/, with given above. The diagonal simplicity constraint is then very easy to deal with. It 

reads: = that is the self-dual and anti-self-dual parts of Bf have equal norm. Thus, there exists a SU(2) 

rotation hf which maps one into the other: 

b_f = -hj' b+f hf (18) 

written in the adjoint representation. Since the bivectors are defined in specific frames, so are the rotations, hfiv) or 
hfij). As explained in [11], the rotation hf has the following interpretation. The Spin(4) group element Hf — [hf, id) 
maps the vector A^^"^ = (1,0,0,0) into a unit vector Nf orthogonal to the Hodge dual of Bf, eijKL ^fBf^ — 0. 
Note that there is a sign ambiguity: we could have simply chosen b-f = hj^b+fhf. Then, the vector Nf = iJ/iV*^"^ 
is orthogonal to Bf. This ambiguity is related to the fact that the simplicity constraints are solved in the continuum 
by two different sectors. The geometric one corresponds to the solution B — ^{e A e) for a tetrad 1-form e, while the 
non-geometric one corresponds to B = e A e. As shown in [ll'|, the sign ambiguity distinguishes between these two 
sectors, although it is irrelevant as far as one only deals with diagonal simplicity. In the following sections, we will be 
able to identify the two sectors at the quantum level, since we will only have to identify them at the classical level in 
the path integral. Our analysis will be done for the geometric sector all along the paper, except when it is interesting 
to compare expressions between the two sectors, in which case it will be explicitly stated. 

Instead of imposing (jlSp with delta functions, we add them to the action with SU(2) Lagrange multpliers g/: 



^diag ^J2^T[Bf{v)Gf{v)j +tr[qfiv){b^f{v) + hfb+fhj\v))\ (19) 
/ 

The choice of g/ G SU(2) is the natural choice since the equations of motion obtained by varying with respect to 
b-f and 5+/ clearly relate g-f and g+z to qf and hfqfhj^. Since the action is linear in the bivectors, integrating them 
projects the fields to these equations of motion, so they will be detailed just below. We allow for a general measure /i 
for the Lagrange multipliers, thus taking into account different possible weights for the quantity {b-f + hj^b^fhf): 

^diag,/,M = J dBf e^"(^/«/) J dhf S,{b-f + h-'b+fhf) (20) 
with J^(a;/) = / dqf fi{qf ) exp {itr[g/a;/]} (21) 

Because of the gauge transformation properties of {b-f + hj^b+fhf), hf transforms with the self-dual SU(2) on the 
left and the anti-self-dual SU(2) on the right, qf transforms under the adjoint representation of the anti-self-dual 
sector, and /i is then taken to be a class function. It is easy to see^ that the strong imposition of the constraints, that 
is 6fj, = 6, is realized for /x(g) = |tr q\. We will find a precise relation between /Lt and the function x of (|17p . 
The integral over the bivector Bf in (pO| can easily be performed with the help of (|A2p : 

^diag,/,M = |tr(g_^^)\r(g_/)| / '^'^f'^'^f ^(^/^ ^(d-fV) ^{d+f Qf hj^) (22) 



* The Haar measure over SU(2) can be related to the Lebesgue measure over TR.^ with the help of a factor [trgl taking into account the 
structure of the multiplication law of SU(2) and its compactness (see appendix . Strictly speaking, it is not true that we have in this 
case S/j, = S since the domain of integration remains compact. One can see however that it has the same effect, by first integrating the 
bivectors. In particular, the trivial weight here gives: maximum value is reached for |6| =0, Ji being the 

Bessel function of the first kind of order 1, up to irrelevant coefficients. 
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The delta functions in this equation are in fact over S0(3) and not SU(2) ^. Thus, the Lagrange multiphers qf stand 
for the holonomies and relate the anti-self-dual part g^f to the self-dual part g^f conjugated by ft./: 

^.a... = J ^i9.f gZ] kj^) (23) 

Since (^5)) imposes g+j and g-f to be related by group conjugation, it is clear that they share the same class angle. 
Since /x is a class function, we have /i(f/_/) = /i((7+/). The key quantity in ([23]), responsible for the non-topological 
character of the model, is obviously S{g^fhfgZj:hJ ), which peaks the amplitude on classical configurations. Notice 
indeed that the group element g-^^fhjgj^ is just the result of parallelly transporting the rotation hj around the dual 
face /. The amplitude thus selects holonomies which preserve the normal vector Nf defined by hf for each face. This 
restriction on the holonomy degrees of freedom can be seen as a characterization or a definition of the model, exactly 
as the BF action serves as an action principle for the model defined by the flatness condition, S{Gf ) 

Using the usual expansion S{g) = 'YljdjXjid) and the Peter- Weyl theorem, the integrations over the rotations 
simply lead to: 



^diagj,,. = H (s+f) ^jf (a-f) (24) 

(trg+/) 



which is precisely PT|) with measure factors related by: /i(g) — (tr (?)^x((7). Thus, to get rid of them in the final 
expression, we have to choose fi{g) — (trg)^, which corresponds to imposing the diagonal constraints, not with a 
delta, but with: 

~ Jl(|f|) + J3(|f|) 

X = 1 < — ^^x) = (25) 

instead, up to irrelevant factors. Here, Ji and J3 are Bessel functions of the first kind defined by: J„(x) = 
IHF 6^^'^°^^ cos{n9)d9. The function (Ji(x) -I- Jz{x))/x is peaked around 0, but with a finite width. This means 
that some fluctuations around the solutions of the diagonal constraints are allowed. These fluctuations are such that 
that integrating over the bivectors Bf yields true delta functions over the group in (j23p . Since we build spin foam 
models from group integrals, this measure is a natural choice. 

More generally, inserting a character ^{g) — Xi{9) the representation of spin I can be translated into the choice: 

4(-^) = ^^^ (26) 

Such a function is not peaked around anymore (it even vanishes at that point for all Z > 1). Instead, it has 
a maximum which is reached for a value of x depending on I. The fact that we nevertheless obtain the BF theory 
restricted to simple representations of Spin(4) can be traced back to the special form of 5^, given as a Fourier transform 
evaluated at + hj^^b^fhj. Then, upon integrating the bivector, the field configurations are always restricted to 

those satisfying g+f = hjg-fhj^. 

Let us introduce a new way to insert the constraints into the discrete path integral. This way of proceeding will 
be argued in this work to be very natural, and will give a very clean action principle for the BC spin foam model. 
Since equation (|22p enforces the SU(2) variables qf to be the holonomies, the constraints should be better introduced 
into the action in a way which respects the group structure and avoids the measure factors (trg_f)'^. Indeed, the 
constraints have been naively added to the action in (fT9|) . so that the latter involves sums of group elements like: 
tT[b-f{g-f + qf)]- The measure factors \tig\ in (|22p are directly due to summing group elements instead of using 
the natural SU(2) multiplication law (see (|A2p ). and do not have any clear meaning. The adapted tool is thus a 
non-commutative addition ® which multiply group elements instead of summing them: 

gtr(&3)©tr(6/i) ^ ^tr{bgh) _^ ^tr{b{g+h)) f^27) 

It deforms the composition of plane waves when the momentum space is given by the group SU(2). This new operation 
is naturally compatible with the SU(2) group multiplication and removes in a natural way all weird measure factors. 



^ Indeed, the expression J d?b exp(jtr(b(;)) does not distinguishes between g and —g (see appendix|A]l. Thus, the expansion into characters 
to be used in the following is: <5so(3)(9) = X]jg]N '^jXjio)^ ^^'^ irreducible representations are of integral spin. 
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More precisely, we will show that it implicitly takes into account the measure (j25p . It is quite convenient when working 
with lattice gauge theories involving elements from both a Lie group and it Lie algebra. In particular, it was already 
shown to arise when looking at observables in the Ponzano-Regge spinfoam model for 3d quantum gravity (soj . 
The diagonal simplicity constraints can be implemented at the level of the action with the following form: 

S.,di,^=Y,TT[Bf{v)Gf{v))®tr[qj{v){b^f{v) + hjH+fhf{v))] (28) 
/ 

= E 9-f{v) Qfiv)) + tr(&+/(f) g+f{v) hf qf hj\v)) (29) 

/ 

Written in this form, the action can be interpreted as a BF theory with 'bivectors'^ Bf oi a, special form: 

Bf{v)^Qfiv)Bf{v), with Qf{v)^{hfqfhj\v),qf{v)) (30) 

This action can be easily seen to reproduce the previous computations, but avoiding the measure factors {tr g^f)^. 
Indeed, the integrations over b-f and only need the formula (|A2p with one of the two group elements being the 
identity. One has instead of ([22]) : 

^*,diag,/.p = / dhf dqf fi{qf) 6{g-fqf) 6{g+f hf qf hj^) 

r (31) 

^Kg-f) I dhf 5{g+fhf g_]hj^). 



■^*,diag,/j — j W^dGr JJ^ ■^*,diag,/,^ with < 



/ 

SO that is recovered with the face weight fc = and the choice fj, — 1. The sum ® can thus be said to naturally 
impose the constraints in a weak manner which is adapted to spin foams, the measure (j25p being taken into account 
through the SU(2) group structure. 

Let us come back to the sign ambiguity of (|18p . Having chosen b^f + hj^b+fhf ~ 0, the field restrictions obtained 
in (j3ip through the BF action are: g+fhfgZ^hJ^ = id. Thus, the sign + between 6+/ and 6_/ has been translated 
into taking the inverse of g-f. This duality gives for the other sector: 



b^f-h-' b+f hf^Q 5(g+f hf g^f h'j^) (32) 



In this case, the orthogonal vectors Nf are left invariant by the group elements ★G/ = (g+/, gZ])- The duality between 
the constraints on the bivectors and the restrictions on the holonomies can be enlightened by using the normals Nf 
instead of the rotations hf. Writing the BF action, for a face /, Tr{BfGf) = ByPf jj, we introduce a Lagrange 
multiplier Xf such that: 



ZdiagJ = JdNfJ dBf J dXf e'ByPf,, + ^Xf (.B,y-'N,.,(*B,),^Nf (33) 

Thanks to the Spin(4) gauge covariance, it is sufficient to analyse this integral for Nf = iV^'^^ . In this case, we have 
= for i,j = 1,2,3, so that only the components i?^' are free. Integrating the BF action then implies that 
pO' — 0, or in covariant form: NfjPj^ ~ 0. Thus: 

^dia..f = I dNf 6(NfjPy^ (34) 



^diagj = j aiy 

which is equivalent to ([?T|) . In contrast, imposing instead NfjBy = leads to the restriction: eijKL Nf Pf^ — 0. 
However, when only dealing with the diagonal constraints, the sign ambiguity is irrelevant, as one can check by 
integrating over hf the r.h.s. of ((32)l . 



The term 'bivectors' is here a misuse of language. Indeed, the following fe^ j do not live in su(2) anymore. We here simply want to 
stress that the action gets a simple form, as a BF theory with bivectors modified, or deformed, by the constraints. 
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III. AN ACTION PRINCIPLE FOR THE BARRETT-CRANE MODEL 



A. Imposing cross-simplicity 

Inserting the solutions Bj = (&+/, —h'^^h-\-fhf) of the diagonal constraints derived in the previous section into the 
cross-simplicity constraints (|12p gives, in vectorial notations: 

b+rb+f'=KfR{hfh-})b+f, (35) 

for pairs of triangles (/, /') in the boundary of a tetrahedron. Note that the sign ambiguity is here irrelevant, as 
expected since the sector of the theory is selected by the choice of a sign when imposing (diagonal) simplicity for all 
bivectors. The constraints ([55)) relate the rotations hf for each pair of adjacent triangles, so that a clear geometric 
interpretation is expected to emerge at the tetrahedron level. Indeed, in [l^, [HI, the constraints are imposed 
into the spin foam model of BF theory by considering that there exists a vector TV,- orthogonal to each tetrahedron 
r (precisely, to the Hodge dual of all bivectors of r). Here, in terms of the rotation variables, the cross-simplicity 
constraints are easily solved by : /i/ h'j} = (jjf where 0/ and (jyf are elements of the U(l) subgroups respectively 

preserving b+f and 6+/'. However, since we are solving (I55t in a given tetrahedron r, the U(l) elements depend a 
priori on r. The one-parameter family which solves the constraints for a given face can be thus be parametrized by: 

b^f = ^hj^ b-^-f hf, with hf = (j)f .^hr and (pf^r b^f \ — b^f (36) 

where hr is taken to be the same for all triangles of r, and is most naturally defined in the frame of r. As the rotations 
hf do not depend on r, there exists a relation between hn and h^-^ when ri and T2 share a face /: 

hrAv) h-^{v)=C^jX{v) <l>f,rM (37) 

This relation holds when hr^ and hT-2 are expressed in the same frame. The geometric meaning is the following. A 
triangle in 4d admits an orthogonal plane, defined by unit vectors N satisfying ejjKL N'^B^^ — 0. On this plane, hf 
stands for the choice of a direction Nf — {hf,id)N^^\ One can check that in turn hr, defined by hr = 4>Jr 
defines a unit vector Nr — {hr,id)N^°'> orthogonal to each face of r, that is eijKL N:}B^^{t) — 0. Notice that these 
(classical) relations correspond to the geometric sector. The fact that hf and hr are related by a U(l) element simply 
expresses that Nf and Nr lie in the same plane orthogonal to /. Similarly, when ri and T2 share the face /, the 
rotations h^ and hr^ are not independent because the normal vectors A^^i and Nr2 are both constrained to lie in a 
well-defined plane, that orthogonal to *Bf. 



A gauge transformation acts on hr as: 



K t> hr = k+r hr kzl (38) 



This specific rule implies that we can always choose a gauge in which hr ~ id, and thus Nr — A^*-°^ being orthogonal 
to all triangles of t. Notice the difference with the rotations hf which only impose simplicity for each single bivector. 
Such a rotation can be gauged fixed to the identity in a given frame, say that of a tetrahedron which has / in its 
boundary. However, it is then impossible to fix the rotations hf oi the other faces of r to the identity. From this point 
of view, the U(l) elements (j)fT enable geometric consistency between gauges of frames and normals to triangles. 

Notice that the relation b^f{T) = —h'Z^b+f{T)hr implicitly contains the elements 4>f,r, as far as a single bivector 
is associated to each face, so that they can be trivially integrated out. This parametrization of the constraints is 
introduced in the path integral measure: 



jX{d%^fSb^f5{\b^f\'-\b^f\') 

xn n s{b+fb+f'^b^fb-f') 

^ (/,/')COt 



J n d%+fd^b^f 1 n n ^(^-/(^) + ^^r' b+fir) k) 



f ^ U,t) 



(39) 

This is in contrast with what can be read in the literature. Indeed, the spin foam models trying to implement the 
cross-simplicity constraints are usually built up only imposing them tetrahedron by tetrahedron. This means that 
the triangulation has been initially broken up into a disjoint union of tetrahedra, so that there are as many copies 
of a given triangle as there are tetrahedra sharing it. A bivector Bfr is associated to each copy of a face /, one per 
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tetrahedron having / in its boundary, and they are all taken to be independent of each other. This corresponds to 
the following measure: 

/n n dH+irdH-fr fUdhr n 5{b.fr + h-H+frhr) (40) 
r fear r ^ r^^ 

We will focus on this specific situation for the rest of the section. However, an essential drawback of this approach 
already appears. The measure P0)l does not encode the relations ([57)1 between adjacent tetrahedra. Indeed, the vector 
Nr defined by hr is only orthogonal to the Hodge dual of the bivectors -B/r, and not to the bivectors labelling the 
other copies of the triangle /. This measure will be supplemented in the following with an action naturally built on 
disjoint tetrahedra and a regluing process using boundary variables for the connection, as in [l^l, leading then to the 
BC model. Since the relations (j37p are not imposed at the quantum level, it will be interesting to know whether or 
not the regluing is sufficiently efficient to make them hold classically. This issue is not clear since we will use boundary 
variables for the connection while the curvature degrees of freedom will not be frozen (unlike in BF theory) . 

It is not surprising that the BC model implements the constraints in this simplified way since an underlying idea 
is the geometric quantization of a single tetrahedron (TtI . [Tsj . As explained above, the same can be said for the new 
vertices. We can also think that restricting attention to the space of intertwiners associated to a single tetrahedron 
in order to solve the constraints, the correlations between tetrahedra may be reduced with respect to ([57)1 . The point 
is naturally that of the gluing process consisting in assigning the same SU(2) representation jf to every copy of each 
face. 



B. A derivation of the BC model: Imposing the constraints on disjoint tetrahedra 

We now propose to study the simplified situation most often considered in the literature, corresponding to ([M)) . 
The aim is obviously to provide a clean link between an action principle and the derived spin foam model, by explicitly 
performing the integrals over the bivectors. It has been argued that strongly imposing the constraints leads to the 
BC model, and that the latter has no equivalent in the other sector of the theory, while the new vertices weakly 
implement the constraints. The main result of this section is that the most naturally obtained spin foam model in 
this framework is the BC model. It in fact corresponds to a weak imposition of the constraints, more precisely to 
the choice of the measure (P5|) allowing specific fluctuations around the solutions of the constraints. Moreover, this 
clean treatment shows that the BC model is indeed a quantization of the geometric sector, and that its equivalent in 
the non-geometric sector is a spin foam model for BF theory restricted to the simple representations of Spin(4), with 
spcciflc face and tetrahedron amplitudes. 

Since we want to impose the constraints independently for each tetrahedron, we assign independent bivectors for 
each triangle of each tetrahedron. This can be achieved by dividing the dual faces into pieces such that each piece, 
denoted (/, r), is completely specified by a dual edge of the face. This means that the triangulation has been initially 
broken up into a disjoint union of tetrahedra. A copy of a triangle is thus specified by a tetrahedron having it in its 
boundary. These variables are very similar to the wedge variables which instead consist of pairs {f,v). Like with the 
wedges, one can build the action for BF theory by coupling the bivectors Bfr to group elements G/r (see figure [5] and 
appendix iBl for the details of the construction). The bivectors BfT are defined in the tetrahedron frames. Important 
ingredients are Spin(4) elements Ljy assigned to the internal edges dividing /. They can be seen as boundary variables 
for the connection, which glue tetrahedra in a certain way discussed in the next section. 

Our starting point is an action for pure BF theory: 

Y,Tv{BfrGfr) (41) 

in which the boundary variables Lfy really reglue the pieces (/, t) of each dual face, expressing the vanishing of the 
curvature. In this broken up triangulation picture, we want to impose the constraints encoded in (j40p : 

b^fr^ -K^ b+fr hr (42) 

which state the existence of independent unit vectors A^,- orthogonal to tetrahedra, or more precisely, ^uklB^^ — 
0. Notice that equation ([1^ is completely classical and corresponds to the geometric sector. As far as integration 
over the bivectors is concerned, we can use the results of section |TT] simply replacing the labels / with /, r except for 
hf which becomes hr- Thus, gathering (|4ip and with Lagrange multipliers qjr G SU(2), we form the action: 

SbI - E Gfr) + Hlfr i^-fr + b+fr K)] (43) 

if-.r) 
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FIG. 2: The left picture represents a dual face, whose edges, denoted r, are dual to tetrahedra and vertices dual to 4-simplices. 
The source and target vertices of r are respectively denoted s(r) and t(r). The dashed lines delimit the pairs (/, r) on which 
we define the bivectors Bfr and which correspond to the triangles of a triangulation initially broken up into tetrahedra. The 
'holonomies' around the pairs (/, r) are defined on the right picture, the base-points being the tetrahedra. The internal links 
carry group elements L , each shared by exactly two tetrahedra. They can be seen as boundary connection variables for each 
T, and integrating them is thought of as a regluing of tetrahedra all together. We have: G/t = Grt{T) ^JtM ^f^M Gs(t)t- 



It is shown in appendix |B] that this action leads to the BC model, using the measure of finite width (P5)) for the 
bivectors -B/t- As explained in section [TTl this measure is related to the fact that group elements are summed in 
the action (|43|) . hke: tr[6_/r(g-/r + 9/r)]. The specific measure (|25|) can be dropped if we use the natural group 
multiplication instead, i.e. the composition © ()27p for each piece (/, t), which implicitly takes it into account. This 
corresponds to the following action: 

^BC = ^^i^-f^ 9~-fT qfr) + tT{b+fr g+fr KqfrK^) (44) 

We have nevertheless noticed that a simpler action leads to the BC model in a very natural way. The underlying 
idea, to say it roughly, is to stick the pieces (/, r) of each dual face using the group multiplication. Indeed, it appears 
in (|44p that the Spin(4) element Qfr defined by Q f^igfr, hr) = {hrqfTh~^ , qfr) is a source of curvature acting at the 
level of the dual edge r. Notice that in the 'time gauge', hr — id, Q^t- is a spatial rotation, whereas we can anticipate 
that it will be a boost in the non-geometric sector (see section UlI Dp . We thus propose to insert these sources into 
the holonomy around each face, starting from the action ^ for pure BF theory: 

S*,BC — ^ Tr^i?/(uo) GygTi QfTi GriviGviT2 Q/ra • • • QfT„ Gt„vo^ (45) 

/ 

where Qfr, naturally defined in the frame of r, is constrained to be of the form: Qfriqfr, hr) = {hrqfrh~^,qfr)- 

S*,BC = ^ tr(^-/(«o) g-voTi qfTi g-Tivi g-viT2 ^/ra • • • ) 

/ (46) 

+tT(b+f{vo) g+vor^ hr,qfr^K! 

a vo,Ti,vi,T2, ■ ■ ■ ,Tn are the ordered boundary simplices of / starting at the vertex of reference vq. The variables 
L fy have thus been eliminated. Notice that in contrast with (|43l) which is well-defined as soon as the faces (and dual 
edges) are oriented, the action l|45p also requires a base-point vq for each face (and thus an ordering). This action a 
priori depends on the choice of the base-points. Indeed, if a bivector is written in another frame using ([3]), the change 
cannot be compensated via as it is the case in pure BF theory, because of the presence of the external sources 
qfr (external here refers to the fact they do not come from a discretization of the connection). For this to work, we 
would need to define rules for parallel transport involving Qfr- However, it will be clear that the partition function 
does not depend on the base-points. 

Due to the hnearity of (|i5l) in the variables Bf{vo), the partition functions projects the group elements coupled to 
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them to the identity: 

Zbc^ [U dGr. [ l[dK n dqfr [ lldBfivo) [] e'T'(^/("o) ^"0^. Q/^i (^/^i ■''^J G.,„, ■.■ ) (47) 

(t.v) r (f^r) f f 

) (48) 

(r.v) T (/,r) / 

It is convenient at this stage to make some changes of variables. First, we redefine the variables g/r to 
g-s{T)r QfT g-Tt{T)7 so that the delta functions on the anti-self-dual sector simply read: 5{(i}ri (lfT2 '■■)• The delta 
functions on the self-dual sector then contains the combinations: hr{s{T)) = gj^s(T)T 9Z-\(t)t ^^"^ ^rit^T)) = 
5+Tt(T) ^'^ 9-Tt(T) which are the result of parallelly transporting hr to the source vertex s(r) of the dual edge r and 
to its target vertex ^(t). Since the group elements g+vT only enters the action this way, we use the Haar measure to 
absorb gZs(r)T ^^"^ ^'^ 9-Tt(T) into g^^j-^-j^ and g+rt(T) respectively. We thus define: 

hT,s(T) = 9+s(t)t hr gZl^r)T 'iT,i(T) = 9+lt{T) 9-Tt(T) (49) 

These SU(2) group elements now define two vectors N^-n^^-^ and N^^s(t)^ both interpreted as orthogonal to r, but seen 
as independent from each other, according to the frame of each 4-simplex sharing it. This change of variables explains 
that from the effective point of view, upon integrating the holonomy degrees of freedom, the BC model deals with 
seemingly isolated 4-simplices. It is equivalent to say that each h-r and 9-tv can be gauge-fixed to the identity. 
The partition function now takes the form: 

) (50) 

= XI / n / n n '^^fdj-f Xi-s (^/ri (?/t2 • • • '?/t„) Xj+/ {hr^,v,, g/n /lri!t,i hr2,v^ ' ' ' g/r„ K^,vo) 

{3±f} {r,v) (/.r) / 

(51) 

In the second line, we have expanded the delta functions over S0(3) into characters, with j+f,j-f G N. Each qfr 
appears twice, so that the othogonality relation (|A3p perfectly works to integrate them. This selects the simple 
representations, j-f = j+f, and recombines into characters the insertions h~\.hri^-,^^vi between the g/rS. To see it, 
notice that the formula (|A3p is simple enough so that we can keep track of the indices. Indeed, the right SU(2) index 
of qfn in the anti-self-dual part is, due to (jA3[) . the same as that of qf^ in the self-dual part, which is contracted 
with /i~^„^/it-2.di • Then, the r.h.s. of the latter is contracted with the l.h.s. of g/T-2 which is also contracted, via (IA3p . 
in the anti-self-dual part, with the right index oiqfr^. This thus forms the character Xjfi^^Tivi^r2.vi)- Let us perform 
explicitly the computation for a triangular face. We are interested in the following quantity: 

/ n ^1- (^1 92 93) X3+ (91 h ) (52) 

3 3 

= E / Dii-J;iq,)Di{-J;{q2)Dg-J;{q,)D^^^^^^^^^^ n^g:!(^-ii/^..+i) (53) 

Oil,2,3 ,/3l,2,3 

= E ]l^'^^,o.A.P^ ^£!(^r/+i hu+i) (54) 

j+ ai,2,3 i 

Q:i,2,35/3l,2,3 

s 

= -^73^ Xj+ (hi2 hi2) Xj+ {^23 ^23) Xi+ (hsi hsi) (55) 
Thus, (|5ip precisely leads us to the integral representation of the lOj-symbol: 

= E n 4 n n ^ n [ / n ^^^-^ n c^^m,. k^^)] (56) 

-EUdlUrr^IlmJf) (57) 
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where u{v) and d{v) are the two tetrahedra sharing / in the 4-siinplex v. 

The action (jlS]) thus provides the BC model with a very natural action of the BF type. The group elements 
/it-,d entering the integral representation of the lOj-symbol represent the unit vectors Nr.v orthogonal to the faces of 
each tetrahedron t when it is seen in the frame of the 4-simplex v. The amplitude then compares the way N^^') is 
mapped into Nt^v for adjacent tetrahedra within each 4-simplcx. Because the change of variables is made on the 
group elements g+vr, the variables h^^v can also be seen as SU(2) holonomies between t and with h-j- and g—^^r 
being gauge-fixed to the identity, that is with flat anti-self-dual transport and orthogonal directions Nr = N^^^^ to 
tetrahedra. 

The obtained amplitude for tetrahedra has already been used in and proposed in [2^ as an alternative to those 
of [1^ and [13] ■ In particular, it does not contain the so-called 'eye diagram'. It corresponds to a trivial gluing between 
4-simplices in the language of [28^1 . Notice however that our derivation uses from the beginning a triangulation and not 
isolated 4-simplices. The decorrelation between the points of view of 4-simplices sharing a tetrahedron comes from the 
integration over the holonomy degrees of freedom. We also emphasize that the method used to impose the constraints 
at the level of the action, for (j44p and (j45p . realizes a natural implementation from the spin foam point of view, which 
can be seen as a weak implementation since it respects the group structure of SU(2) entering the BF action. It is 
shown in appendix IB] that, from the point of view of the action (j43p . it corresponds to introducing the constraints 
into the path integral with the measure (PS)) . which is not surprising in the view of section |TT1 As for insertions of 
observables, it is clear that the actions (|44p and produce the same results for insertions of holonomies. As for 
B-dependent observables, one has to be careful since both actions are defined with different bivectors. 



C. Equations of motion for the BC action 

Having shown that the BC model can be derived by imposing the constraints on disjoint tetrahedra, one can still 
clarify the situation described by the model, and in particular the gluing process, by looking at the equations of 
motion of the discrete action. Instead of (|45p which requires an explicit ordering of the dual edges, let us consider 
the action (|44)) which treats tetrahedra more symmetrically. We obviously have: 

q/T = 5l/r and g+fr K gZ]^ K'^ = id (58) 

which still hold at the quantum level since is linear in the bivectors. These can be thought of as characterizing 
the model, exactly as 5{Gf ) defines the BF spin foam model and 6{g+fhfgZ^hJ^) restricts to simple bivectors. If we 
interpret G/r as parallel transport around / seen in r, the right equation of (j58p means that N^- is left invariant by 
parallel transport. In the time gauge, it says that Gfr is simply a spatial rotation. Notice that in the continuum, 
the equation of motion w.r.t. the field B still involves B (since the constraints are quadratic in B). However, the 
analogous equation is here independent of the bivectors. 

To extremize the action w.r.t. group elements such as qjr, we use right invariant vector fields given by Sqfr qj^ G 
su(2). Extremizing the action w.r.t. the multipliers g/r (and using (j58p ) naturally gives the simplicity constraints 
(|42p . The equations of motion for Grt(r) and Gg^T^^T are respectively: 

^ Bfr = and ^S- ^f- = (59) 

/car /car 

These relations are naturally consistent with (|42)l and ((58)) . In agreement with the remarks below ([49]) , these relations 
for the self-dual sector only, together with (|42p and imply those for the anti-self-dual sector. The first equation 
of ([55)1 is the so-called closure relation and can be seen as a discretization of the continuous equation (IaB = 0, 
integrated over a tetrahedron of the triangulation. For pure BF theory, the second equation of ([55)1 is equivalent to 
the closure relation since Gfr = id. However, due to the sources g/r, this is not the case anymore and it means that 
the closure relation is maintained when 'parallelly transporting' each of the four bivectors around the corresponding 
triangle with Gfr- 

Equations ([i^ and ([55)1 ensure metricity for each tetrahedron independently. Let us thus study the regluing based 
on the variables Lfy. The standard spin foam model for BF theory is often built using some division of the dual faces 
into plaquettes p, which allows the assignment of independent variables, like bivectors, to each piece p of every face. 
Such a division corresponds to breaking up the original triangulation into a disjoint union of simplices. For instance, 
in the specific situation here studied, the triangulation has been initially broken up into a disjoint union of tetrahedra, 
so that a copy of a triangle is specifie d by a tetrahedron having it in its boundary. The wedge variables are similar, 
but constructed with pairs (/, w) (see [20] for instance). Then, we can assign independent bivectors Bfr to the pairs 
p = (/, t), which are related by group variables, here the elements Lfy. These latter depend on the dual face and live 
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along the boundary shared by a pair of plaquettes belonging to the same face. We emphasize that they should be 
seen as a tool for gluing simplices rather than a tool for parallel transport between frames. Such a gluing is however 
a non-trivial operation. A priori, the two following actions for BF theory are different: 

J2t^t[Bp{v)Gp{v)) ^J2^T(^Bf{v)Gf{v)) (60) 
p f 

The fact that they both lead to the same spin foam model is due to the specificity of BF theory, that is the triviality 
of the holonomies. For adjacent plaquettes of a face /, the bivectors are related by: 

Bf,r2 = G^/^^ Bf^n GriT2 (61) 

when the action (|4ip is extremized w.r.t. to Lfy. Such a relation is necessary for consistency between frames and is a 
discretization of the continuous e.o.m. cIaB = which states that B is constant up to parallel transport. Moreover, 
implementing it all along the boundary of /, it implies that the holonomy G/ leaves Bj invariant, Gj Bj GJ^ = Bf. 
The fact that the degrees of freedom of G/ lie in 11(1)6^^: x \J{l)b_j, is a basic property of Regge calculus. 

We expect this relation and (pT|) to be modified in the presence of a source for parallel transport. Indeed, the 
variables q/r sue responsible for a new relation when extremizing w.r.t. Lfyi 

%,T2 = Q7,T2 '^ri\-2 Bf,Ti Gr,T2 Q f ,T2 wlth Q fr = {hrqfrK^^lfr) (62) 

This is precisely the relation which makes the action (|l5)) independent of the base-points. This specific relation 
induces a non-trivial effect on the relation ((37)) between the normal vectors of neighbouring tetrahedra. Using simplicity 
under the form (H^ together with ([5^ . one obtains that hr2{T\)h~^ {t\) = g+riT2^T2g-\T2^T2 leaves 6+/ri invariant, 
as expected (or equivalently, the Spin(4) element HT^r2iji) = (''■r2 ('''i)^ti"'^('''i)i ^T/('''i)'ir2 ('''i)) preserves the bivector 
Bfr-^). But transporting this relation in the frame of T2 with (|62p shows that (t2)/i^^^(t2) = ^T23-Tir2'*r/ff+Ti7"2 
does not leave fe+/.r2 invariant, due to insertions of qf^. We have instead: 

hT2 qjr2 &+/,r2 {hT2 qjr2 K^ (^2))^^ = &+/,r2 (63) 

with b+fr2 = hr2qfT2K2^ b+fr2 hr2qJr2K2 (64) 

or equivalently, (/i^^ qj^^ K^i'^^.), qjr2 hr2{T2)) acting on Q/^^ B/^^ Q/r2 gi'^ss B/^^- 

To conclude this study, it is clear that the action (weakly) imposes cross-simplicity on disjoint tetrahedra. 
The variables Lf^ have to reglue them together, but precisely because of simplicity, we do not obtain the expected 
relations when parallelly transporting the bivectors describing the same triangle. The variables g/v play indeed the 
role of holonomies but with a dependence on the face. This in turn alters the correlations between neighbouring 
tetrahedra. The rotation /i^^/i"^ is supposed to preserve &+/, since the two normals live in the same plane orthogonal 
to However, this relation is not satisfied in every frame. 



D. In the non-geometric sector 

The explicit implementation of the classical relation (|42p into the path integral shows that the BC model is a 
quantization on isolated tetrahedra of the geometric sector of Plebanski's theory. The non-geometric sector corresponds 
to: 

h^fr = ""^ b+fr (65) 

Since — tr(6g) = tr(feg^^), it is sufficient to replace each qf^ with its inverse in the self-dual (or equivalently the 
anti-self-dual) part of the action to obtain this sector. Working again on disjoint tetrahedra, we propose: 

S =^ tr{b-fr g-fr qfr) + tr(6+/^ g+fr K qj^ h:;^) (66) 

However, the model defined in this way turns out to be ill-defined. It indeed cruelly depends on the parity of the 
number of boundary edges of each dual face. This difficulty can be overcome by imposing cross-simplicity not simply 
on isolated tetrahedra, but instead for each tetrahedron in each 4-simplex separately. The details of the construction 
are given in appendix [Cl A tetrahedron r being shared by exactly two 4-simplices, s(t) and t{T), it means that 
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t{r) 



Gfr 



frv' 



f,r,t{r)) 
if,r,s{r)) 



FIG. 3: It represents a dual face, whose edges, denoted r, are dual to tetrahedra and vertices dual to 4-simplices. The dotted 
lines delimit the wedges, that are pairs (/, v) and which correspond to the triangles of a triangulation initially broken up into 
4-simplices. The dashed lines represent the breaking up of each 4-simplex into disjoint tetrahedra and form the half-wedges, 
consisting in triplets (/, r, u). The standard spin foam model for BF theory can be defined with these half- wedges, with the 
help of holonomies Gtv and of 'gluing' variables L/r and L/„ living on the internal links. 



the existence of the unit normal Nr is separately imposed in s(t) and t(T), with the same rotation hr, but with 
different multipliers qj .^ g(^^-^ and qj: .,. f(^^y The flipping of q in the self-dual part of the action can be recorded in 

Q frv — {hr iJtv Qfrv)- This Spin(4) element is a boost in the 'time gauge', while it was a spatial rotation in the 
geometric sector. 

5"= ^ TT(^Bfrv GfrvQfrv^ (67) 

(f,r,v) 

The assignment of the variables is illustrated in figure [31 The plaquettes (/, r, v) will be called half-wedges, since 
they arise from superimposing the wedge division and that into pairs (/, t). Since a tetrahedron is shared by two 
4-simplices, the number of copies of a triangle is twice that obtained with the pieces (/, r). Thus, with regards to the 
above-mentioned parity problem of the number of dual edges, this construction selects the case where each face has 
an even number of edges, since all happens as if tetrahedra were chopped into two pieces. Notice that this framework 
could have also be used in the geometric sector, leading again to the BC model. In particular, the e.o.m. relating the 
variables Bfrv are of the type (|62p and the variables Qfrv can be naturally composed to form: Qfr — Q fTs{T)Q fTt(T)i 
which is consistent since the normal to r is encoded in the same rotation h-r for s(t) and ^(t). 

Like in the geometric sector, we can 'stick' tetrahedra thinking of the Q frvS as insertions of curvature source into 
the holonomies around the faces at each half- wedge, and replace each Q/r in (|45p with Q^t-^j-^) Qj^^^y. 

S* = Tr^i?/(uo) Gygn QfTivo Qfrivi GtiVi G-vj^r2 Qfr2Vi ' ' (68) 

/ 

This gives, in self-dual/anti-self-dual decomposition: 

Si, = ^ y ^^f^-/ ff— "01"! QfTiVo QfTiVi 9 — TxVi 9 — VlT2 QfT2Vi 

f ' (69) 

+ .9+t)ori 'Z/rii'o ^7ti"i 9+viT2 K2 <7/t2«i ' ' ' ) 

as an equivalent of (|45p in the non-geometric sector. In the geometric sector, qj^, is simply replaced with qfrv in the 
self-dual part, so that they can be obviously composed to form the multiplier qfr = qfTs{T)<}fTt{T)- However, this is 
not possible anymore in the non-geometric sector because of the flipping of the multipliers in the self-dual part of the 
action. 

The steps of the computation are similar to those of the derivation of the BC model. First integrate the bivectors: 

-UqTi QfTiVo QfTiVi 9-TiVi 9-ViT2 9/T2Ui ' ' ' j 

(t,v) T {f,r,v) f (70) 
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and expand the resulting expression into characters. Then integrate over the variables qf,T,v, thus selecting the simple 
representations. To see more precisely what happens, let us keep track of SU(2) indices while using (|A3p . The right 
index of g-vori is contracted with the left index of qfnvo- Thanks to (jA3p . the latter is equal to the left index of 
qjrivi self-dual sector. Using (|A3[) again, this index equals the right index of Qfrivi in the anti-self-dual sector, 

which is contracted with g-nvn and so on . . .This process forms the holonomy g-f{vo). In the self-dual sector, the 
rotations hr disappear and we also end up with g^f{vo). The remaining integrations are taken to be over the group 
elements Gr = Gs(T)TGrt{T) allowing for parallel transport between 4-simplices, with Gf = Grj^Gr2 ■ ■ ■ ■ Taking into 
account the dimensional factors dj^. , we have: 

[l[dGrJ2Il rr^—f - n 4 x., (5+/) x,, (5-/) (71) 

= E E 114 H n \^ ni5jSpin(4)((j/,J/);(Hr,^-.)) (72) 

The model can thus be straightforwardly interpreted as a spin foam model for BF theory restricted to the simple 
representations of Spin(4). However, in comparison with the model obtained in section |TT] with the diagonal simplicity 
constraints, the face and tetrahedron amplitudes are different and are the direct effects of the constraints^ ([55]) . 



E. About the Immirzi parameter 



The Immirzi parameter 7 represents a one-parameter ambiguity in classical pure gravity. It can be included into 
the usual Palatini-Cartan action by using the Hodge dual -k of so (4): add to the curvature F the quantity -kF with 
coefficient 7: 

^Hoist = J eijKL Ae^A (i^^^ + 2j{kFf^) (76) 

This is the so-called Hoist action [sij. The equation of motion for the new term in this action is simply: c?yi(eAe) — 0, 
which clearly does not modify the original equations. Notice that it is equivalent to leave F unchanged and proceed 
toeAe— >eAe + 7*(eAe) instead. In Plebanski's formulation, the field B is then similarly transformed, but the 
simplicity constraints also need to be changed. We thus prefer to use the version (|76p with B = ★(e A e) together with 
the usual simplicity constraints. 

Let us define for g = cos 6* id 4- isinOn ■ a E SU(2) and a e M, the SU(2) element whose class angle is a9 and 
direction n: 

:= cos (ad) id -I- i sin (aO) n ■ a (77) 
For integral a, this is obviously the group multiplication. The BF action with an Immirzi term reads: 

S.y= I (l + 7)tr(6+ AF+) + (1 - 7) tr (5_ A F_) (78) 



Let us briefly present a model based on SU(2) BF theory in which face or tetrahedron amplitudes are affected by adding some information. 
Let <f>f and xfif be in SU(2), and add their commutator <l>fipf(pJ^^J^ to the usual discrete action: 

S = J2^^{bf9f<l>fi^f'i>j'i^j') (73) 
/ 

The interpretation of the model is as follows. The e.o.m. are: 

gf ipfipf <f> J '^ip J = id ^j:bfij>J^=^fbf^f = bf (74) 

together with the closure relation. Thus <j)f and ipf give a source of curvature which is classically constrained to live in the U(l) subgroup 
preserving bf. In this case, <p-f and tpf commute, so that Qf = id. However, at the quantum level, g/ is not constrained to be the 
identity. Although the partition function does not project on flat connection anymore, the vertex of the resulting spin foam model is 
still the 15j-symbol: 

Z=T,U^U^^^(^f) (75) 
O/} / - 

Only the face weight, instead of dj , takes into account the modification of the action. In this regard, notice that it is well-known 
(since the works of Ponzano and Regge) that the face weight dj^ is required for the model to be triangulation independent. 
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and is clearly equivalent to pure BF theory, but it will not be the case anymore in the presence of the simplicity 
constraints. The issue is now to introduce 7 at the discrete level. Consider the su(2) element g±f standing for the 
discrete curvature, that is the holonomy of the connection along the boundary of the dual face /. In a given chart, 
consider that the edges of / are of order e. When e goes to zero, we have: g±f w 1 + e^Fj-^j, where -F±|/ is the 
component of F± along the face /. Thus, defining 7± = 1 ± 7, it is clear that: « 1 + e^7±-F±|/, and a natural 
choice to reproduce the continuous action (|78p is: 

S.^Y. i^+f Sl^f) + {b-f (79) 
/ 

Let us study the imposition of the diagonal simplicity constraints into this framework, to compare the result with 
the usual ansatz for the introduction of the Immirzi parameter into spin foams [H, Il5| . Diagonal simplicity still 
expresses the fact that the bivectors are simple, and is solved by relating 6+/ to &_/ through a rotation hf which 
represent the choice of an vector Nf orthogonal to *-B/. At the classical level, the Immirzi parameter modifies the 
relation between the self-dual and anti-self-dual parts of the curvature. We have indeed: 

b^f^^hj'b+fhf ^ gl+^ = hf gl-f hj' (80) 

which means that the vector Nf is not preserved anymore by the Spin(4) element Gf = (.g+/,g_/), as it is the case 
in ((23|) . but by the element (g!j.p(?Ij). Notice that the quantities j± are not relevant in themselves, but rather their 

ratio is. Indeed, ((80)) is equivalent to: g+f = hf g'^Sf'^^ /ij^, saying that g+f lies in the conjugacy class of g^^j'^^ ■ 
We can then easily read the limits 7 — s- ±00: 

5+/ = hf g-_) hf (81) 

which indeed corresponds to the non-geometric sector. For a general 7, since the action and the constraint are linear 
in the bivectors, ([50)) still holds, avoiding some measure factors like in section HI) 

/ n ^Gt. Y^\dhff> (5+/ hf gl-/''^ hf) (82) 

{t,v) f 

= E / n dGt. n (.9+/) {ai-/'^ ) (83) 

The resulting spin foam model is however ill-defined for a generic 7: g-f is a product of group elements g-tv and 
we do not know how to apply the operation ([77]) on each subfactor. We thus ask for a quantization of 7 such that 
rij = 7-/7+ = (1 — 7)/(l + 7) is an integer. The usual ansatz for introducing the Immirzi parameter is based on the 
canonical analysis of the constraints, in which the bivectors are quantized as the generators of so (4). This leads to 
the following sum: 

E / n ^G*- n ^».^/ (5-/) (84) 

fa/} it,v) f 

with the same quantization of 7. If the dimension dj = 2j + 1, instead of the spin, is scaled by in ()84p . the 
difference with ([55)) can be seen to be a measure ambiguity, and both are easily related. Indeed: 

Xk,,Jg) = = X(n-i)/2(5) X]i9 ) (85) 

for kj^n = nj + (n ~ l)/2 chosen such that: dk ^ — ndj. We thus propose to include the Immirzi parameter with the 
following action: 

tT{b+f g+f hf qf hf) + tr(6_/ gT'f Qf) (86) 

/ 

where qf G SU(2) is a Lagrange multiplier. Then, we can think of ([5^ starting from this action with the following 
weight for qf: X(n-,-i)/2('?^^"^)7 determined by the Immirzi parameter. To deal with cross-simplicity, the basic idea 
is thus to replace gZf in ()44[) with g_f ■ The computations are quite cumbersome (in particular, the number of 
appearance of l-fv depends on 7) and the result is not enlightening, so we do not reproduce them. 
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FIG. 4: The left picture depicts the sequence of couphngs for a given face. The Unks carrying the measure representations 
kf^T are the boundary edges of /, dual to tetrahedra. It is a very natural generalisation of the BC model in which all fc/,^ 
are taken to be zero. The right picture represents the vertex of the model. The faces (the links in the picture) are labelled 
with irreducible representations of Spin(4). These are intertwined at each tetrahedron (node in the picture) with 

the representations kf^r which are shared by others vertices. 



IV. SOME GENERALISATIONS OF THE BC MODEL 



We present in this section two types of generalisations of the BC model. The first type is a natural extension due to 
choosing a non trivial measure for the Lagrange multipliers qfr- It does not change the physical setting of the model, 
but gives a family of models in a simple way. Furthermore, the vertex amplitude of these models naturally exhibits 
the fusion coefficients of the new spin foam models [l3, [HI, ■ The second generalisation consists in weakening the 
imposition of the constraints. The new spin foam models can be interpreted in this framework. 



A. The fusion coefficients 



We have seen that the action (|45p implements the constraints naturally with regards to the SU(2) structure and 
gives a straightforward derivation of the BC model, with a trivial weight for the variables qfr- Let us still consider 
the action (|45p , which imposes the existence of normal vectors independently for each tetrahedron, and introduce 
a typical gauge invariant weight ^J^iqfr) — Xkfrilfr), where the spins {kfr) are initially given. Omitting the integrals 
over the rotations h^- and over the holonomies G^t, the partition function is: 

^{kf.A = / n '^1frXkf,^{qfT)Y[s(^g^r^r^qfT^9-T^T2■■■qfTS) ^ (/lr„^5+r„Ti ^-n ^/n /iri^g+nra ^r^ ' ' ' g/r„ ) (87) 
if,r) f 

if vo,Ti,vi,T2, ■ ■ ■ ,Tn IS the sequence of boundary simplices around /. The element of parallel transport Gnn+i is 
naturally Gnn+i — GnviGy^n+i ■ Expanding the delta functions into characters assigns representations j+/ and j_/ to 
each face. The elements qjT now appears thrice. Thus, we no longer have j+/ = j-f, but instead these representations 
are intertwined at each tetrahedron with fc/r according to (|A5[) . 



in which the rotations have been absorbed on the right of the group elements g+vr (or equivalently gauge-fixed to 
the identity), r and r' denote the two tetrahedra sharing / in w. The sequence of couplings for a face is described in 
figure m 
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The typical vertex is easily derived. The integrals over the group elements g±Tv give rise to 4-valent intertwiners 
*±T,t), between the four representations j±f labelling the faces of r. The SU(2) indices carried by v are contracted 
to form a 15j-symbol for Spin(4): 15jgpin(4)((j+/, {i+rv,i-Tv))- On the other hand, the tetrahedron amplitude, 
depicted in figured is due to contracting the indices living at each t. The links carrying the representations fc/r, 
coming from the intertwining of j^f and J-/, relate the 4-simplices v and v' which share t. Remarkably, this result 
can be re-expressed in terms of the fusion coefficients /-^ j_ used in the spin foam model of [l^ . Indeed, one only 
needs to insert the identity on the space of intertwiners jj^ (g) jj^ (g) j/3 (g) jf^ ^ (D as a sum over a complete orthogonal 
basis consisting of 4-valent intertwiners labelled by the internal representation 1^-, as in figurelH The partition function 
then reads: 

with 

Ay ^ E 15jSpin(4)((j+/, j-/); (j+r^,«-T.u)) (j+/,j-/,A:/r) (90) 

These models thus form a very natural generalisation of the BC model, involving the same basic idea to implement 
the constraints. At first sight, these models seem to be close to the unflippcd model of 'l^. The differences are the 
following. One should first impose j^f = j+f- Second, having obtain a vertex similar to that of [12], the main point 
is how the amplitudes Z^fc^^i, should be summed. The physical meaning of the representations kfT- is unclear and a 
general measure is a complete expansion into characters Xk;^ ■ For instance, the vertex amplitude is weighted in 
[T^ l with the Clebsch-Gordan coefficients CjJ_^-^g^^. However, such a coupling between the face representations and 
those coming from the expansion of the measure is here not accessible. 

The amplitudes ((90)) obtained in this section show that the fusion coefficients fl^ ^_ naturally arise. Indeed, the 
simplicity constraints explicitly relate the variables of the self-dual sector to those of the anti-self-dual sector, the 
Spin(4) covariance being preserved by hr- This is achieved in the functional integral through the Lagrange multipliers 
qfr- From this point of view, the BC model does not show up these coefficients because it is the simplest model in 
this class, defined by the weight 1 for each qfr- When choosing all fc/r to be zero, the 3jm-symbols of ([88|l become: 

—5j_^j^5a.-on SO that the intertwiners i±Tv can thus be summed to form the well-known BC 
intertwiner ibc- 



B. Weakening the constraints 

A second proposition to look at the BC model with a more general p ersp ective is to weaken the imposition of the 
constraints. This is indeed the main idea from which the new vertices of |lOl.[ll| have originated, and which have been 
proposed to remedy the drawbacks of the BC model. The main approach makes use of the Hamiltonian analysis, first 
to quantize the bivectors, then to argue that the cross-simplicity constraints, seen as second class constraints in this 
framework, have to be weakly imposed. 

We have already argued that the physical framework of the new vertices is that of the BC model. We thus still 
impose the constraints tetrahedron by tetrahedron, and search for implementing the constraints in a weaker sense. For 
this purpose, we use, instead of the variables qfr, different multipliers for the self-dual and the anti-self-dual parts of 
the action (|45p . q+fr and q-fr, which are coupled through gauge invariant functions p{q-f-fr, q^fr)- This corresponds 
to the following partition function: 

f Y[ dq+frdq-fr p{q+fr,q-fr) fY[dBf e»'K''-/("«)9--o-l9-/..i9-l..l-)+»tr(''+/(-o)s+„o-l''-l9+/..l'»;">+xi^l-) 

(91) 

= Y\. '^l+f-rdq-fT p{q+fr,q-fT)Y\_^(^9-voTiq-f,Tig-TiVi- ■ S(^g+y„rihr-,q+f (92) 
if^r) f 

again omitting the integrals over hr and G^r- The gauge invariance of p reads p{krq+frkr^ , krq-frkr^) — 
p{q+fr,q-fr) for each pair (/, r) and kr £ SU(2). The previous class of models is recovered by choosing a cou- 
pling of the special form: p[qj^-,q^) = 5{qj\-qZ )pi{q+) for the geometric sector. Notice however that one can now 
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switch between the geometric and non-geometric sectors only with a change of measure p{q+, 9-) = p{<l+, 9-^)- This 
is indeed equivalent to replacing q-fr in the action with it inverse. Thus, the boundary between the two sectors is 
weakened by the weakening of the constraints. 

Since p can be expanded into matrix elements of g+ and q_ contracted with intertwiners, let us restrict attention 

to p{q+fr,q-fT) = D^All'l^B+f-^i'l+f-r)D^All'l'B-;S^-f-r)- O^c Can then perform the integrals over q±fT using the 
orthogonality relation (|A3p . which enforces the equalities j+/r = J+/ and j-fr — j-f- Upon absorbing the rotations 
hr into the holonomies g+vT as in the previous models, we have: 

(t.v) f T i-i-fCdT 3 + f 3- J f^Q^ 

where r and r' are the two tetrahedra sharing the triangle / in the 4-simplex v. It is thus clear that the function p 
directly encodes the couplings of the model. The features which do not depend on the choice of p are the labelling of 
faces with Spin(4) representations {j+f,j-f), the face amplitude dj^jdj_^, and the fact that the vertex is a product 
over the triangles of / of functions of g^l^g+vr' and gZlrg-vr' ■ 

Let us now show that ((93|) is a general setting allowing to understand the new vertices as weakly imposing the 
constraints. Indeed, we are already assured to be able to write the vertex with the fusion coefficients. Moreover, in 
the representations labelling faces are now exactly those coming from the expansion of p. We thus look for the 



above-mentioned coefficients C^J . Consider the special choice: 



= J (i'^k+b^iXi^^-k-b^) (94) 

3 

where \j,n) are SU(2) coherent states, here associated to pairs {f,T). They have been introduced in [2^ and further 
used because their semi-classical behaviour eases the geometrical interpretation of the constraints. However, they do 
not seem to naturally appear in our computations. The function ([M)) . which leads to the flipped model of ^] (or the 
unflipped model by choosing p{q^,qZ^)), can be compared with that corresponding to the BC model: 

PBc{q+,q~) ^6{q+qZ^) = ^d] j d^h {j,n\q+ qZ^\j,n) (95) 

3 

which shows that ([M)) indeed weakens the natural SU(2) measure corresponding to the BC model. 
Expanding (j94[) into matrix elements, we obtain the following expression: 



{3f}-^ U,r) f r . 



(96) 



with ^'f = y ]X '^G^^ ]X^^f''^f-^'\9+r'v g+vr\jf,nfr) {jf,nfr'\g-T'v g-VT\jf,nfr) (97) 
rCdv f 

This expression has been shown in [l^ to reproduce the flipped model of The other sector, that is the unflipped 
model of [HI, can be obtained with a simple change of measure which leads to: 

^«=y" n '^'^"^ Yl{jf,nfr'\g+r'v g+vr\jf,nfT) {jf,nfr\gZlr gZl>Jjf,nfr') (98) 

rCdv f 

An important part of the content of these models has been here introduced by the measure p, used to weaken the 
constraints. The other physical inputs, that is the choice of the variables, and the imposition of the constraints on 
disjoint tetrahedra, are common features shared with the BC model. From this point of view, a better understanding 
may require a precise analysis of the function p in (j94p . Using {j,n\g\j,n) = {cos0 + isinO u ■ n)'^^ for g ~ cos0 + 
i sin9 u ■ a, the integral over the unit vector n can be explicitly performed with repeated use of the binomial formula. 
p{q+, q-) then appears as a function of the class angles of the SU(2) group elements q+ and q_, and of the dot product 
of their directions. The resulting expression seems however quite cumbersome. 



V. PHYSICS OF THE BC MODEL AND BEYOND 



It is the first time, as far as we know, that the BC model is derived using a path integral formulation on a lattice. 
We thus focus in this section on the physical setting of the model. The variables have clearly identified geometric 
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roles. The clean treatment of the situation points out the drawbacks of the BC model with regards to quantum 
gravity, which can be traced back to the issue of gluing between tetrahedra, and thus to that of their correlations. 

The BC model has been originated by considering an isolated 4-simplex, and imposing the constraints on each 
of its tetrahedra. It has then been thought that the failure of the BC model to describe quantum gravity resides in 
this specific fact, and that imposing the constraints on a complete triangulation leads to the new vertices, depending 
on the sector, [12,]. We have however seen here that the BC model can be thought of on a triangulation (discarding 
boundary terms), without isolating 4-simplices from each other, and that the new spin foam models can also be 
considered in this framework. A rotation h-r is associated to each tetrahedron and geometrically represents the 
direction Nt = (ft,T-,id) t> 7V*^°) perpendicular to t, without refering to the 4-simplices sharing r. The fact that the 
BC model deals with manifestly decorrelated vectors TVt-.-u arises from the change of variables (|49l) . The rotations 
are defined at each tetrahedron, but the amplitude, comparing the ways two adjacent tetrahedra map iV*-"^ into their 
normal directions, takes place at each 4-simplex. Thus, integrating over the holonomy degrees of freedom decorrelates 
hr{s{T)) and hr(t{T)). This is in contrast with the derivation of the BC model proposed in 0, [12]) in which the 
constraints are imposed 4-siniplex by 4-simplex (and then tetrahedron by tetrahedron). Moreover, the form of the 
simplicity constraints used in |12| to derive the new vertex is the same as that which has been here shown to reproduce 
the BC model when explicitly integrating bivectors. 

It is clear by looking at the equations of motion of the BC action that the BC model contains some correlations 
between neighbouring tetrahedra ri and T2- Such correlations are given by ([63| on-shell, while the kernel of the BC 
model (I56p is due to comparisons between the rotations ft-n and h^^- However, we have seen in section Hill that the 
latter are directly related by the simplicity constraints. Indeed, the cross-simplicity constraints are solved by assigning 
rotations h-r to tetrahedra, but without forgetting that for tetrahedra sharing a triangle, these normals are correlated 
since they live on the plane orthogonal to Therefore, correlations between tetrahedra are intimately related to 

the gluing process. These correlations are not correctly taken into account in the BC model, because of the gluing 
and of the fact that the rotations hr are considered as independent variables. 

The correlations are written in (j37p in terms of elements 4>f,T living in the U(l) subgroup preserving the self-dual 
part &+/ of each bivector. The U(l) subgroup preserving 6+/ is of major physical importance because it enables 
to relate the normal directions of neighbouring tetrahedra, and thus describes the whole set of orthogonal vectors 
to a given triangle. Let us introduce bivectors BfT like in section IIIIl Then we can choose a normal A^/,r of 
reference for each tetrahedron sharing /, encoded in hjr G SU(2), to impose diagonal simpHcity for each piece (/, r): 
h^fr = ~^Jt ^+fT ^fT- Cross-simplicity is very easy to impose in this framework: 

hfr = e"'f-^^+f--ihr (99) 

with a/r S [0,47r). 6+/r is the direction of b+fr, so that hfrh^^ preserves 6+/r. We then need to introduce some 
explicit rules for parallelly transporting Bfr between adjacent tetrahedra. The norm |6+/r| = |^-/,t| is the same 
for all T, and Grr' sends the directions {h+fT,h-fr) of Bfr to those of S/t'- One can show that, as far as only the 
directions are concerned, it is equivalent to give such rules for the self-dual parts 6+ fr together with a rule to relate 
the vectors Nf^ between tetrahedra sharing /. Not surprisingly, the latter is of the form ((57|) : 

hfr h-jl,{T) = hfr 9-rr' hj^, g'^^, = e^^' ^+/-^ (100) 

with i^rr-i G [0,47r). Again, the r.h.s. leaves &+/t invariant. It means that the issue of parallelly transporting Bfr is 
intimately related to that of consistently imposing simplicity. 

Another problem of the BC model, as already mentioned in the study of the equations of motion of the BC action, 
is that it does not reproduce any discrete analog of Einstein's equations. We expect in particular the following 
relation: G f{T)BfrGj^{T) — Bfr, where Gf is the holonomy around the dual face /, as a discretization of the 

fact that {e^ A e"') A (*F)jk = in the continuum (or from the e.o.m. (IaB^'' = as already discussed). This 
relation between the triangles and the holonomies around them can be seen in terms of the vectors N fr- Indeed, 
Gf{T) G V{l)b+f^ X U(l)b_j,^ implies that though Gf{T) does not leave Nfr invariant as in (|58p. it sends it instead 
onto a new vector which still lies on the plane orthogonal to *i?/(r): 

eijKL [G/(r) > Nfr] '' B^^^ = (101) 

The analog of Einstein's equations in the Plebanski formulation is obtained by extremizing the action w.r.t. the field 
B. In the BC action, this gives equations (|58p which are independent of bivectors and require that the normal Nr 
is left invariant by Gfr- This is due to the fact that the action is built on disjoint tetrahedra and is linear in the 
bivectors. However, the above discussion shows that the directions and the norms of h±fr have quite different roles. 
A study of a discrete action built on these ideas and of the resulting spin foam model will appear in a coming paper 

il. 
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Conclusion 



We have shown that the Barrett-Crane spinfoam model can be derived from a discrete action principle and is 
related to a quantization of the geometric (gravitational) sector of Plebanski's theory. It is based as usual on a BF 
term discretized on a triangulation. The simplicity constraints are independently imposed on each tetrahedron and 
we explicitly perform the integrals over the bivectors. The Barrett-Crane spin foam amplitude is recovered for a 
specific choice of measure weakly imposing diagonal and cross-simplicity. This measure can be translated into the use 
of a non-commutative composition of plane waves, and it is argued to be the natural way to proceed in spin foam 
computations. We have also proposed an action which sticks tetrahedra together using the SU(2) group structure. 
This method straightforwardly leads to the same results and thus implicitly takes into account the above-mentioned 
measure. 

It is known that the BC model does not precisely correspond to gravity since it does not reproduce the expected 
correlations in semi-classical computations of graviton type [9| . The Englc-Pereira-RovcUi (EPR) and Freidel-Krasnov 
(FK) models have been proposed to remedy the drawbacks of the BC model. We have here obtained a Lagrangian 
formulation of these models by weakening the implementation of the constraints in the lattice path integral. In our 
view a better understanding of these models requires a careful study of the regluing process, which sticks tetrahedra 
together after imposition of the constraints. The regluing which enables to recover the BC model is that usually 
considered in pure BF theory and unconstrained BF-like theories [20j . However, the Lagrange multipliers can be 
seen as a source of curvature, so that this regluing yield incorrect rules for parallel transport of bivectors. This 
issue is related to the geometric consistency of the constraints and to a faithful accounting for correlations between 
neighbouring tetrahedra. These correlations are detailed in section |Vl Equation ([57]) , whose trace gives the dot 
product of the normals to adjacent tetrahedra and thus contains the dihedral angles of Regge calculus, describes them 
in terms of variables living in U(l) x U(l) subgroups of Spin(4) which leave each bivector invariant. A proposition 
using precisely such variables is under preparation p32| . 

In the non-geometric sector, the amplitude is built on the 15j-symbol for simple unitary irreducible representations 
of Spin(4), exactly like the model only implementing diagonal simplicity. However, they both differ due to the 
amplitudes associated to triangles and tetrahedra. A precise introduction of the Immirzi parameter into the partition 
function is also proposed and shown to differ from its usual introduction in a measure factor. 

Finally, the fusion coefficients, at the root of the EPR and FK spin foam models, are derived starting from the 
setting of the BC model supplemented with a non-trivial measure for the Lagrange multipliers. Thus, they appear 
to be very natural objects in the context of spin foam quantization of gravity. From this point of view, the main 
questions concern the spin data labelling the links of the boundary spin networks and the summation of these 4-simplex 
amplitudes, that is to say the choice of weights for triangles and tetrahedra. 
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APPENDIX A: SOME USEFUL FORMULAS 



SU(2) group elements are parametrized hy : g = cos9 + isinO n ■ a, where a are the Pauli matrices satisfying 
[ai, aj] — Cfc. We define the projection of g onto the Pauli matrices p — sm9 h. Thus, the (unnormalized) Haar 
measure over S0(3) (restricting the parametrisation to 9 E [0, ^]) is related to the Lebesgue measure on : 

dgso(3) = n ^ (Al) 

The factor (1 — jP')^ restores the compactness of the group when working with the p variables, \p\ < 1. 



Using test- functions, it is easy to show that : (5so(3) (ff^ ^) = y 1 ^ Pg ^'''^\pg ~ Ph)- Rewriting the delta over R'^ 
with the usual Fourier transform, we have : 

^ Pa 
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where b — Vri is a Lie algebra su(2) element, integrated over with the Lebesgue measure, with = — fcTj- Such 
a formula does not admit a straightforward generalisation as soon as more than two group elements are involved. 
Indeed, a sum of p variables, phi + Ph^^ is not in general the projection of a group element onto the Pauli matrices. 
Equation ()A2|) is very simple, but given the parametrisation of the constraints introduced in section it is the only 
formula needed in this work to explicitly perform the integrals over the bivectors. 

After the integrals over the bivectors, one encounters integrals of products of matrix elements of SU(2) representa- 
tions. The basic formulas are that expressing the orthogonality of these matrix elements, 

? D^ab\9) D'f^^^ig) = -^6,,,, 5ac5M. (A3) 

and that giving the decomposition of a tensor product of representations, 

= E Ca'^A Ci^^i D^Alia)- (A4) 

J,A,B 

with \ji — j2\ < J ^ ji + j2 and —J< ji,j2 < J, and where the coefficients C^u^ J g^j.^ Clebsch-Gordan coefficients 
between the irreducible representations ji (g) j2 — > J- One can easily show, using (|A3|) and (|A4p . that: 



J dg D'^^^\g) D^^{g) D'^f*{g) = C^T^ (A5) 

j dg D[i^\g) Dl^\g) D^^* {g) D^^^ ig) = E ^^a^?^ (^^"^'b C'^'.'a C^H (A6) 



J,A,B 

One can express these relations with the help of the 3jm-symbols which enjoy better symmetries than the Clebsch- 
Gordan coefficients thanks to: 

mim2m3 \ i v ji \jn^ jn2 —Tils J 



APPENDIX B: ANOTHER DERIVATION OF THE BC MODEL 



We show how cross-simplicity can be imposed on disjoint tetrahedra in the most naive way, that is by simply adding 
them to the action. In this framework, the BC model corresponds to a specifically weak imposition of the contraints 
corresponding to the measure ([25]), and not (|40)l . This feature is naturally hidden when we stick the BF action with 
the constraints using the SU(2) multiplication, as in section ITTTl a process which respects the SU(2) structure of the 
spin foam. 

To impose the constraints independently on each tetrahedron, the triangulation has been first broken up into a 
disjoint union of tetrahedra. Then, we associate bivectors Bfr to each face of each tetrahedron, and consider them to 
be all independent from each other. We will then impose the cross-simplicity constraints by asking for the equality of 
the rotations associated to the triangles of a tetrahedron. 

Let us first build the BF action using the variables Bfr- The construction is similar to that using the wedges. From 
the point of view of the 2-complex dual to the triangulation, a wedge is a pair (/, v) formed by a face and one of its 
vertices (i.e. a triangle seen in a particular 4-simplex, in the triangulation language). The variables (/, t) that we use 
here are pairs consisting of a face and one of its boundary edges (i.e. a triangle seen in a particular tetrahedron): draw 
segments from the center of each dual face towards each of its vertices, as shown figure [2] We thus use decorrelated 
bivectors Bj^ for each pair {f,T). The BF action then couples them to some 'holonomies' Gfr around the pieces 
(/, r). Two pieces (/, ri) and (/, T2), for tetrahedra ri and T2 sharing / in the same 4-simplex v, have a common link 
delimiting them (see figure [2), and which carries a Spin(4) element, Lf^y entering the elements Gfr^^- Integrating 
these variables in the path integral can be seen as a regluing process. As discussed in section IIII CI this regluing 
perfectly works as far as there are no sources of curvature, which is the case in BF theory but not in our action for 
the BC model. The following action reproduces the usual spin foam model for BF theory: 

SbF = E Tr(i?/r Gfr) (Bl) 

if.r) 

As noticed in section IIIIl this action is a priori different from the standard discrete BF action ^ . The particularity 
of BF theory that it is stable under breaking up the triangulation into tetrahedra, 4-simplices or both is directly due 
to the fact that it deals with trivial parallel transport around dual faces. 
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Let us now present the computations of the partition function from the action Sbc^ in which the constraints for 
the geometric sector are naively introduced: 

Sbc - J2 ^/O + Tr[g/,(6_/> + b+fr K)] (B2) 

if.r) 

where -B/t, Qf,T and hr are defined in the local frame of r. Sbc is a function of the Spin(4) group elements Gy-r 
(collectively denoting the basic variables Gs(t)t E^nd Grt{T)) and Lfy, the SU(2) Lagrange multipliers qj^, the rotations 
hr standing for the normals to each tetrahedron, and of the bivectors -B/r- The partition function thus reads: 

Z^= f l[ dGyr n ^^/- / n^'^- / n ^^/- e'^"" n S^^ib^fr + KH+frhr) (B3) 
iv,T) if,v) T (/,r) (/,r) 

with / dg/^ ^(g/r) exp {i tr[g/^ x/r] } (B4) 

JSU(2) 



As tetrahedra have been unstuck, integrating over the bivectors is very easy with (jA2[) and makes the measure 
factors |tr((7)| appear. Then, the integrals over the Lagrange multipliers relate the self-dual parts of the elements 
Gfr = {gj^fr, S-fr) to thcir anti-self-dual parts: 

Z^= [ U dGyr n dLfy [lldhr n 7^7^^ '^(3+/- 9Z]r K') (B5) 

Note that , similarly to (|23p , the group element 17+/^- ft,,- 5 ~ is a kind of parallelly transported GfrOhr of the rotation 
hr, so that Gfr has to preserve the vector Nr- This restriction on the fields is responsible for the non-topological 
character of the model, and is independent of the measure fi of the Lagrange multpliers. It could have been in fact 
anticipated from the action since the latter is linear in the bivectors, so that equation (jB5|l simply expresses the 
projection onto the corresponding equations of motion. 

We now restrict attention to the special choice ([^5]) for the measure fi. This amounts to considering only the delta 
fimctions over the group in (jBSp . and equivalently to use the specific function Sf^sci^) — '^^'^^^^'^^^^'^^^^'^ for each pair 

(/, t), instead of the strong condition 5^^\x). Integrating the boundary connection variables Lfy is then easily done 
since each of them only appears twice. In particular it implies the equality of the representations j/r, assigned to 
the pairs (/, r) via the expansion of the delta functions into characters, to a single representation j/ for each face. 
A triangle in a single 4-simplex 1; being shared by exactly two tetrahedra, let us denote them u{v) and d{v). The 
partition function is thus: 

zbc= [u dGry / n E n 4 n rr^—r n [ n x.. {Ki.^ («) h-^i^ ^] im 

The rotation hriv) stands for the parallelly transported of hr to the 4-simplex v along the boundary of a dual 
face. In (|B6[) . hr is only transported to its source vertex (in the dual picture) s(t) and its target vertex ^(t): 
hr{s{T)) — gj^s(^r)T hr gZl(^r)T ^^'^ hr{t{T)) — g^^^^^-^ hr g-rt(T)- It mcaus that the vector Nr, orthogonal to the 
tetrahedron r and defined in its rest-frame, is parallelly transported to the two frames of the 4-simplices sharing 
it. Then the character Xjj{hu{v){'^)h~^(y){v)) compares the way the vector iV^"^ is mapped into Nr for two adjacent 
tetrahedra at the 4-simplex v. Since the self-dual components g+vr of the group elements realizing these transports 
only appear here in the amplitude, g-vr and hr can be absorbed into them by the change of variables (I49p . This 
directly leads to the integral representation of the lOj-symbol, (|5^ 



APPENDIX C: THE BC MODEL IN THE NON-GEOMETRIC SECTOR 



The model turns out to be ill-defined if we simply think of imposing 6-/t — h^^b+frhr = by changing qfr with 
its inverse in the self-dual part of the action (|45p. It indeed crucially depends on the parity of the number of dual 
edges for each dual face. The model becomes well defined if we use the wedge variables, superimposed to the variables 
(/, r) previously used, according to figure[31 Let us stress that all the models derived in this work can be reformulated 
using these half-wedge variables. However, since they are only necessary in the present case, we have chosen not to 
use them before to simplify the constructions and to obtain more straitghforward interpretations. 
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Let us break up the triangulation into a disjoint union of 4-siniplices. Thus, there exist as many copies of a given 
triangle as it is shared by different 4-siniphces. In the dual picture, it is represented by a division of each dual face into 
pieces corresponding to the choice of a boundary vertex, that is into pieces (/, v) called wedges. Let us now break up 
each 4-simplex into tetrahedra. Since a tetrahedron is shared by two 4-simplices, it means that the number of copies 
of a triangle is twice that obtained with the wedges. Thus, with regards to the above-mentioned parity problem of 
the number of dual edges, this construction will select the case where each face has an even number of edges, since all 
happens as if tetrahedra were chopped into two pieces. A copy of a face is specified by the data of a tetrahedron and 
a 4-simplex containing it. The corresponding division of a dual face is presented in figure [31 and we call each piece 
(/, T, v) a half-wedge. 

We are now looking for imposing the constraints on each half-wedge, but with the same rotation /i,- for half-wedges 
of the same face and same tetrahedron, that is imposing the existence of a unique normal N^- such that N^-jB^j'^^ = 
for the two copies of r. In the naive way to impose the constraints on the action, bivectors have to be independently 
assigned to each copy of a triangle, -B/rw, and the regluing of the simplices can be partially done with boundary 
variables associated to the links delimiting the pieces {f,T,v) (the dashed and dotted lines in figure [3]). But this can 
be done more efficiently with the arguments leading to ([^5]) . Having multipliers qfrv to impose the constraints on 
each half-wedge, these can be viewed as sources Qfrv — {hT<lfTvh~^ ,qfTv) inserted into the holonomies Gf for the 
geometric sector. The flipping of the sign from (fT8|) to |32|) translates into the change qfrv qj^v self-dual 
part of the action. We thus define Qjtv = [hrq'jl^h~^ ,qfrv)- Obviously, for the geometric sector, one can define 
the variables qfr = q^^^j^^gy^jj^) and Q/v = Q fTs{T)Q frtir) appearing in the action of the BC model (|45p . This is 
however impossible in the non-geometric sector, because of the flipping qfrv —>■ QJtv self-dual part of the action. 

Let us mention a subtlety of the construction. Since a tetrahedron exists in two copies in the shattered triangulation, 
these copies have a priori different frames: one for r in the 4-simplex s(t) and one for t in ^(t). We thus introduce 

a Spin(4) element Gr''^''*'^'' responsible for the parallel transport between these two frames. The rotations describing 
the same normal in the two frames are related by: /i^''^^'''' = g^^^*'"^'' h'^'"^'''' g'^^^*''^^~^ . This relation supplements the 
following action: 

SC-S = ^^"^(j^-f 9-VoTi qf,Ti,Vo 9-Ti <lf,Tl,Vi g-TiVi g-V^T2 qf,T2,Vi ' ' 

/ (CI) 
+ tr(&+/ g+vor. h^T^qjX^^^ g^^ qj^^^.h^",^^ ' g+rr.n hil^^jX^^^---) 

This action being linear in the bivectors, integrating them in the partition function is easy: 

(C2) 

Expanding the delta functions introduces representations j+/ and j_/ labelling faces. Integrating the SU(2) vari- 
ables q/Tv is very standard, using the orthogonality relation (jA3p . It enforces the equalities = j_/ = j/, eliminates 

(v) 

the rotations hr , and recombines the holonomies around the dual faces Gf = GvgriG^°^^ Grj^y^ • • • , exactly as ex- 
plained in section UlIDI We can thus turn to the variables Gr = 6*5(1-) G'r'"^^*^^''GT-t(T-) representing parallel transport 

between the 4-simplices s(t) and ^(t). It means in particular that the elements Gt^^^**'^'' can be gauge fixed to the 
identity, so that the action can be simplified to that presented in section lllli ([55]) . One finally ends up with the 
spin foam model of equation ([7^ . whose difference with the model for BF theory restricted to Spin(4) simple rep- 
resentations, implementing simplicity of each bivector independently and presented in section |TT1 lies in the face and 
tetrahedron amplitudes. 
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